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Abstract 

We obtain global space-time weighted-L 2 (Morawetz) and L 4 (Strichartz) estimates for a 
massless chargeless spherically symmetric scalar field propagating on a super-extremal (over- 
charged) Reissner-Nordstrom background. To do this we first discuss the well-posedness 
of the Cauchy problem for scalar fields on non-globally hyperbolic manifolds, review the 
role played by the Friedrichs extension, and go over the construction of the function spaces 
involved. We then show how to transform this problem to one about the wave equation on 
the Minkowski space with a singular potential, and prove that the potential thus obtained 
satisfies the various conditions needed in order for the estimates to hold. 

1 Introduction 

There have been many studies of the well-posedness and decay of scalar fields in a given space- 
time whose metric satisfies the Einstein equations of general relativity, both as a precursor 
to the study of the stability of that metric, and as a means of probing various censorship 
conjectures. Although our main goal in this paper is to obtain estimates for scalar fields, we 
need to address the issue of well-posedness as well. 

1.1 Well-posedness of the Cauchy problem for scalar fields 

This question arises in studying the phenomenon of singularities in general relativity. The 
classical notion of a singularity of space-time is that of geodesic incompleteness. The weak 
cosmic censorship conjecture states that generically, singularities of spacetime must be hidden 
inside black holes, instead of being "naked", i.e. visible to distant observers. The strong form 
of this conjecture posits that generically, spacetimes must be globally hyperbolic, i.e. possess 
a complete spacelike hypersurface such that every causal curve in the manifold intersects it at 
exactly one point. Global hyperbolicity ensures that the spacetime has deterministic dynamics, 
since there is a Cauchy surface whose domain of dependence is the entire spacetime. In the 
case where the space-time is not globally hyperbolic, which can happen for example if the 
metric is singular 1 on a time-like curve, well-posedness of the Cauchy problem for the scalar 
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x The minimum regularity required of the metric for the local existence and uniqueness of geodesies to hold is 
C 1,c * , a > 0, so that any point in the manifold where the metric fails to be C 1 '" needs to be removed, which can 
easily result in the loss of global hyperbolicity We note however, that the singularities of metrics considered in 
this paper are in fact much stronger, and the curve-integrability conditions of |5j do not hold for these examples. 



wave equation has been suggested [2U as a substitute for geodesic completeness in determining 
how singular the spacetime actually is, and perhaps to see if quantum effects have any chance 
of regularizing the dynamics and restoring predictability |1(J| . This is particularly relevant if 
censorship conjectures somehow fail to be true and naked singularities turn out to be more 
abundant than otherwise allowed by them. 

Let us now recall the set-up from 24 : ,: Let (Ai,g) be a static, stably causal space-time, i.e., 
one that admits a hypersurface-orthogonal Killing vector field whose orbits are complete 
and everywhere timelike. A massless scalar field on Ai satisfies the equation 

g^VuV^ = 0. (l) 

Suppose we specify initial data for i/i on a hypersurface S that is everywhere orthogonal to T^. 
If Ai is not globally hyperbolic, £ will not be a Cauchy surface and data on £ will determine 
ij) only on the domain of dependence D(Y,). The aim of |24| was to define a physically sensible 
recipe for determining tp everywhere in Ai. To this end, one rewrites Q in the form 

dfip = aD a (aD a ip) 

Here a = y/—g(T, T), t is the Killing parameter, and D a is the covariant derivative of the 
Riemannian metric induced on S. One then views 

A := -aD a (aD a ) 

as an operator on the Hilbert space 

H = L 2 (S,a _1 do-), 

where da denotes the induced volume form of S. The point of this definition is that with 
respect to the inner product of TC, A will now be symmetric and positive, although not yet 
self- adjoint since the initial domain of A would have to consist of sufficiently smooth functions. 
If we take this initial domain to be C£°(£) (smooth functions of compact support on S) then 
A is also densely defined, and the equation we want to solve is 

(d 2 t + A)^ = 0, V| E = /, = 9- (2) 

It is then a classical result (20], that the above problem is well-posed, provided one replaces 
A in this equation with (one of) its self-adjoint extension(s) Ae- Such self-adjoint extensions 
are guaranteed to exist for real, symmetric operators [2Hl- The extension may not be unique 
though, which would imply that there is still an ambiguity about the dynamics. This may be 
interpreted as having to specify boundary conditions for the scalar field "on the singularity". 
One finds that there are three proposals for removing the ambiguity and restoring determinism 
to the dynamics of scalar fields: 

(1) One can of course restrict attention only to the cases where the self-adjoint extension 
is unique (the so-called essentially self-adjoint case), and declare that spacetimes where the 
operator A is not essentially self-adjoint are "quantum-mechanically singular" ^01 (meaning 
the singularity remains even if quantum particle dynamics, i.e waves, are considered in place of 
classical particle dynamics, i.e. geodesies). In this approach, the naked singularities present at 
the center in both the negative mass Schwarzschild (m < 0) and the super-extremal Reissner- 
Nordstrom (|e| > m) spacetimes are quantum-mechanically singular, and nothing more can be 
said about the evolution of scalar fields on these backgrounds. 
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(2) Another possibility is to characterize the singularity that leads to non-unique self- 
adjoint extensions for A as having a U(N) "hair," jllj where N is the common value of the 
two deficiency indices n± := dim(ker(^4* ± i)) of A. This is because on the one hand there is 
a one-to-one correspondence between the self-adjoint extensions of A and unitary maps from 
ker(^4* +i) onto ker(A* —i) [23], and on the other hand, as shown in any "reasonable" way 
of defining the dynamics on the whole spacetime, would necessarily have to arise from some 
self-adjoint extension of A. We will show that the naked singularity of the super-extremal 
Reissner- Nordstrom solution is in this sense, quite hairy (see Remark 0). 

(3) A third approach is to distinguish one self-adjoint extension from among all the pos- 
sible ones, as being somehow more "natural" or "physical". It was suggested in ]23J that the 
Friedrichs extension of A, i.e. the one coming from extending the corresponding quadratic 
form, is such a natural choice. We note that the Friedrichs extension of A always exists, and is 
unique by construction 2 . Moreover, as was shown in |2j, the Friedrichs extension is the only 
self- adjoint extension of A under which the resulting dynamics of (J2J) agrees with that of the 
corresponding first-order (Hamiltonian) formulation, 



obtained via the (unique) skew-adjoint extension of the corresponding operator h. 

In this paper we will take the third approach, and work with the Friedrichs extension of 
A. This is mainly because finiteness of the H 1 norm is used several times in the course of the 
proof of our estimates, and among the extensions of A, the Friedrichs extension is the only one 
whose domain is contained in H . It may be that some of our estimates hold for the other 
extensions of A as well, but our proof does not extend to those cases. 

1.2 Estimates for scalar fields 

The next natural question to consider for scalar waves after well-posedness is obtaining esti- 
mates for them. This is relevant among other things to the question of stability of the metric 
as a solution of the field equations, since expanding the perturbation functions in tensor har- 
monics yields a sequence of scalar wave equations |17j , and having good estimates for each one 
of them seems necessary-although not sufficient-for proving linear stability of the metric under 
a small perturbation of its data. Moreover, the question of stability has an obvious connection 
to cosmic censorship if the solution to be perturbed happens to have a naked singularity. 

Most of the results obtained so far in the direction of estimates concern either boundedness 
of the field ^3 [23 ^ or its dispersion, i.e. the decay rate, with respect to the foliation 
parameter, of the supremum of the field on the slices of a given time-like or null foliation of 
a portion of the space-time ^3 [21 El El HI HI El El E] In this paper we take a different 
approach to this problem, by obtaining estimates that bound the global space-time L -norm 
(and weighted-L 2 norm) of the field, in terms of appropriate norms of its data on a given Cauchy 
hypersurface (see 013 for a similar approach). Such Strichartz (resp. Morawetz) estimates, 

2 The claimed non-uniqueness of the _ff 1 -based extension (which is equivalent to the Friedrichs extension) of 
this operator in the case of superextremal Reissner-Nordstrom, 11 §IV.A.2], is due to the authors' error in not 
realizing the Sobolev space H 1 as the closure of functions under the corresponding norm (see Section 3). 
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as they are known in the literature, are a staple of the modern theory of partial differential 
equations, and have proved very useful in the analysis of nonlinear evolution problems. We 
will be looking at the simplest non-trivial situation possible, namely, a spherically symmetric 
scalar field on a given spherically symmetric and static background space-time 3 . Furthermore, 
the background here is allowed to have a naked singularity, the main example we have in mind 
being the super-extremal (naked) Reissner-Nordstrom solution of the Einstein-Maxwell system. 
Our main result is 

THEOREM 1 Let (M,g) be the Reissner-Nordstrom space-time manifold, with mass m and 
charge e, such that 

|e| > 2m, 

and let tp be a massless, chargeless, spherically symmetric, scalar field on M . Then there exists 
a constant C > (which may depend on e and m) such that 

W\ L ^ g ) + Wp'^Wl^) < CE 1/2 (V), (3) 

where p denotes the area-radius coordinate on Ai, defined in Section 2 and E 1/ / 2 (V') *s the 
conserved ^-energy of the field, defined in Section^ 

The outline of this paper is as follows: In Section 2 we introduce the family of static spherically- 
symmetric Lorentzian manifolds, define a global system of coordinates on them, and write down 
the evolution equation satisfied by a scalar field on such a manifold. In Section 3 we introduce 
the notion of self-adjoint extension that is necessary to make that evolution problem well-posed, 
and define the function spaces that are going to be used. In Section 4 we state the various 
conditions that need to be satisfied by the metric coefficients of a manifold in the family we are 
considering, such that a scalar field on it will satisfy the estimate (J3J). This is then proved by 
transforming the problem to one about the flat wave equation with a potential, appealing to 
the result in [5], and transforming back to the original problem. The last section contains the 
proof of the fact that the metric coefficients of the super-extremal Reissner-Nordstrom solution 
satisfy the above mentioned conditions. This is accomplished by reformulating the problem in 
the language of real algebraic curves and appealing to the compactness result of |22j . 

2 Scalar fields on static spherically symmetric 
backgrounds 

2.1 The space-time 

Consider a four-dimensional connected spherically symmetric static space-time {M,g). More 
precisely on M we assume a time-like action of M and a space- like action of SO (3, M) commuting 
with it. These actions should be without fixed points, except that at most one M-orbit is allowed 
to be SO(3, M)-fixed, in which case it will be called the time axis. We restrict our attention to 
the cases where the SO(3, M)-orbits off the time axis are spheres. 

3 It is possible to extend our result to fields that are a finite sum of spherical harmonics, as in but we do 
not pursue it here since this approach is not likely to yield an estimate for the general, nonsymmetric case. 
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2.2 The coordinate system 

Let T = T{p) denote the Killing vector generating the E-action at p € M. and let A = A(p) 
denote the area of the 50(3, E)-orbit passing through the point p G M. We define two 
E x SO(3, ]R)-invariant functions on M as follows: 

a := y/-g(T,T), p := ^Aj^. 

The quotient space (Q,g) ■= {M, g)/ 50(3, E) is a two-dimensional Lorentzian manifold. 4 We 
can assign coordinates t, r on Q in such a way that 

g ab dy a dy b = a\r) (-dt 2 + dr 2 ) . 

The Killing field is clearly just dt in these coordinates. We pull back t and r to M. It is easy 
to see that we can assign angular coordinates Q = (9, <f>) on M. in such a way that the metric 
becomes 

9llv dx» dx v = a 2 {r) {-dt 2 + dr 2 ) + /?(r) 2 (^ 2 + sin 2 (j)d6 2 ), (4) 

and the action of SO(3, M) is the usual action on the unit sphere. These coordinates are unique 
up to translations in the t, r coordinates and rotations in the 6, cf) coordinates. If there is a time 
axis we can then arrange that r = there. 

We note that in these coordinates the volume form of the spacetime is 

d^ g = a 2 p 2 dCldrdt, 

while the induced volume form on the Cauchy hypersurfaces t = const, is 

da = ap 2 dQ.dr. 

Here 

d£l = sin (pdOdcf) 

denotes the volume form on the standard 2-sphere. We also note that there are at least two 
other volume forms on these hypersurfaces that are natural to consider, namely 

/ 2 2 // 1 2 

da = ada = a p dSldr and da = oT da = p dQdr. 

The significance of da' is that L p -norms based on it behave as one would expect with respect 
to the t-foliation: For / : M — ► E and all p > 1, 

ll/l|LP(d/x s ) = II Il/Ilx^(dcr') ||x^(dt)- 

The significance of da" will become clear in the next section. 
4 Possibly a manifold with boundary if there is a time axis. 
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2.3 The wave equation 

Consider the action functional A = \ j M Ldp g corresponding to a real massless scalar field 
ifi : M — ► M. The Lagrangian density is 

L = g^V^V^. 

In the above-given coordinates 

L = -a~ 2 i/j 2 + oT 2 ^ 2 . + p~ 2 \ J7^\ 2 , 
where JJ denotes the unit-spherical gradient, and 

|^| 2 = l^| 2 + -4-|5 e <A| 2 . 

sin <p 

A scalar field ip that is a stationary point of the action A, subject to a given set of initial values 
(ipo, -01 ) on the Cauchy hypersurface t = 0, satisfies the following Cauchy problem 

dfy + AiP = Q, 0(0) = ^o, ^0(0) =Vi, (5) 

where 

A := -^d r (p 2 d r ) - ^4L, (6) 
p z p z 

with ft. denoting the Laplace-Beltrami operator on the unit 2-sphere. Note that the operator A 
is symmetric and positive definite with respect to the inner product given by the volume form 
da", i.e. for 0,-0 € C c °°(S t ), 

/ 0AV> da" = [ ^Acp da", 

JT, t 7E t 

and 

v2 



/ 0,40 d<r" = / 02 + ^| X70| 2 da" > 0, 



with equality only if = 0. 



3 Self-Adjoint Extensions, Energy, and Sobolev Norms 

As explained in Section ^ the evolution problem (JSJ) is only meaningful if the operator A 
there is replaced by (one of) its self-adjoint extension(s) Ae- In the event that the self-adjoint 
extension is not unique, the one we choose to pick is the Friedrichs extension Ap, obtained by 
extending the quadratic form naturally associated with the operator A, namely 

Q A (0):= / r 2 + 4l^0| 2 d<r". (7) 

Thus the Cauchy problem we are actually studying is the following 

d 2 ^ + A F ^ = 0, V>| E =/, d t ^ s =g. (8) 
Let T^y be the energy tensor of the field, defined as 

Tfj, u = V u ip - ^g^L. 
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In particular, in the above given coordinates 

Let X be a time-like Killing vector field for M. Then the one-form P defined by P(Y) ■= 
T(X, Y) is divergence-free, i.e. * d * P = and thus *P is a conserved current. In particular, 
let X = dt, then Pq = Too an d Pi = Tqi, so that the only nonzero component of *P is 
(*P)i23 = p 2 Too, and the conserved quantity is the energy 

E[$ := / p 2 T 00 drdft = f T 00 da". 

Thus using the dc" volume form, Too is identified with the energy density of the field. 

It therefore makes sense to also use da" to define Sobolev spaces on the E^'s: Let H 1 (E t ) 
denote the completion of smooth compactly supported functions on Ej with respect to the 
norm 

11/11^:= / \fr\ 2 + ^\Wf\ 2 da". 
We define H (E^) in the same way, i.e. completion with respect to the norm 

ll/l&o := / |/| 2 da". 

The Sobolev spaces H s (Et) for < s < 1 are then defined via interpolation between the above 
two spaces, and one uses duality to define them also for — 1 < s < 0. We note that by these 
definitions, 

and moreover 

m = \ (ll*l( Et ) + WW*)) ■ 
More generally, for < s < 1 we can define the s-energies 

E,M :=5(WH.(E t ) + H^llH.-i(E t )) J 
and it is not hard to see that they are all conserved under the flow @: 

4 Transferring to Minkowski Space 

It is clear that if we can view (jSJ) as an evolution in Minkowski space, we may then be able to 
apply known theorems in order to get the estimates we want. This is of course only possible 
if the manifold M has at least the same topology as M 4 (or M 4 minus a line, in order to 
allow a singular time axis). We will be making this assumption from now on, and denote the 
coordinates on R 4 by the same letters as those on Ai, namely (t,r, O) G R x R + x § 2 , and a 
t-slice in R 4 is still denoted by Ej. Let us define u : R 4 — > R by 

u{t,r,fl) := wX/>(t,r,Sl). 
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Note that 



u 



\u\ 2 r 2 drdVL 



\ 2 da" 



2 

H°(St) ' 



(9) 



'St JT.t 

We can easily check that for ip S C£°(Et) we can perform an integration by parts to obtain 



u 2 + Wu\ 2 + 4^r(-)]V + -^(-)^(^) 
pi r z p r p 



r 2 drdn 



a 



r 2 drdfl =: Qb(u), 



where the "potential" V is defined by 

V(r)=p"(r)/p(r). 

Let B denote the operator whose associated quadratic form is Qb defined above, i.e. 



(10) 



We then have that 



c) + A = -{d 2 +BY-. 

p r 



Note that the first term in the definition of B coincides with the radial flat Laplacian. In 
particular, if ifi satisfying © is spherically symmetric, i.e. ip = ifj{t,r), then u can be viewed 
as a radial solution of the flat wave equation with a potential 



uu-u Tr u r + V(r)u = 0, u(0, r)=/(r), u t (0,r) = g(r), 



(12) 



with / = pipo/ r an d g = pipi/r. 

Since by (|1U|) the quadratic forms corresponding to A and B are equivalent, we can identify 
the Friedrichs extensions Ap and Bp of A and B, as well as the Sobolev spaces based on 
them, i.e. let Ti. s denote the S-based Sobolev space, defined to be the completion of smooth 
compactly supported functions on M 3 \ {0} with respect to the norm 



W 



■■= \\(B F y/ 2 f\\ 



L 2 ■ 



We then have that for Isl < 1 



\u\\W 



This is because the 1? spaces agree by © while the H l spaces agree because of QlOjl. On the 
other hand, if we let 

P := -A + V = -^d r (r 2 d r ) -\fi + V, 

then it is clear that on the subspace of radial functions, Bp coincides with Pp (the Friedrichs 
extension of P), and so do the corresponding Sobolev spaces, which we will denote by TC ra ,d- In 
particular if we define the s-energy of u to be 



S s [u] := \\u(t)\\ H s + |K 



^rad 
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then 

£ s [u]=E s [^} (13) 

and it is conserved by the flow (|12j) . 
In jHj it was shown that solutions to 

(d 2 + P F )u = (14) 

satisfy certain weighted-L 2 (Morawetz) and IP (Strichartz) spacetime estimates given below, 
provided the potential V meets certain criteria. In the case of a radial potential V = V(r) on 
M 3 these criteria reduce to the following three conditions on V: 

supr 2 V(r) < co, (15) 

r>0 

inf r 2 V(r) > -1/4, (16) 

r>0 

sup r 2 ^-(rV{r)) < 1/4. (17) 
r>0 dr 

The version of the main result in [5] for radial potentials (and general data) is as follows: 

THEOREM 2 Let V G C 1 ((0, oo),R) satisfy Q3Q^7^, and Zet P := -A + V(M) w/iere A 
is i/te Laplace operator with domain C£°(M. 3 \ {0}), and let Pf denote its Friedrichs extension. 
Then there exists a constant C , depending only on the quantities on the left in U5llb]ll\ ), such 
that any solution u of satisfies 



" 1 «||l2(r4) + ||u||l4 (R 4) < C£y 2 [u]. (18) 
Now, the space-time LP norms of ip and u are related in the following way: 

p-2 

a 2 [- 

I M J-ooJS 2 J0 



P 

LP(M) 



2 -) \u\ p r 2 drdndt. 



In particular 



Similarly 



1*<m> = / (y) M-Adt. (19) 



We thus have proved the following theorem regarding estimates for the scalar field if) on Ai: 

THEOREM 3 Let Ai be a Lorentzian manifold that is homeomorphic to R , admitting a 
timelike R action and a spacelike SO(3, R) action commuting with it, in such a way that exactly 
one WL-orbit, called T, is SO(3, W)-fixed. Let (i, r, O) 6lx R + x S 2 6e i/ie coordinate system on 
Ai as in Section^ with T = {r = 0}. Let g be a Lorentzian metric on Ai that is of the form 
0), is C 3 outside T, and is such that the functions p and a satisfy the following conditions 

(i) sup r>0 (r 2 y) < oo 

(ii) inf r>0 (r 2 F) > -1/4 



9 



(iii) S u Vr>0 (r 2 £(rV)) < 1/4 

(iv) inf r>0 (-^) > 

ar 

where 

V ( r > := ~T2- 
p dr z 

Then there exists a constant C > 0, depending only on the quantities on the left in the conditions 
above, such that any spherically symmetric solution of 

dfy + Apip = 

satisfies 

Wp^WlHM) + UWiHM) < CE 1/2 [# (21) 
Proof: We simply observe that by 



1 1 

\\P~ */>\\i? + H^llr* < -^V^uWtf + T/^IHI^ 4 

where d denotes the quantity on the left in condition (iv). The result then follows from Q18j) 
and JT3J). 



5 Super-extremal Reissner-Nordstrom 

The Reissner-Nordstrom manifold (M,g) is the static, spherically symmetric solution of Einstein- 
Maxwell equations. It is characterized by two parameters: mass m and charge e. It can be 
shown that for a metric of the form (jlj to satisfy the Einstein-Maxwell system, one must have 

and 



/ 2m e 2 

a = \l + -=. 23 

V p p 

It thus follows that when |e| < m the function p is bounded away from zero and there is no 
time axis. This is called the sub-extremal (black hole) case. When |e| > m the manifold A4 
has the same topology as R 4 minus a line, the time axis is at r = and the metric g is highly 
singular there. This is the super-extremal (naked) case of Reissner-Nordstrom. 

The ODE 1)22(1 can be solved to express the "tortoise" coordinate r as a function of p. In 
the super-extremal case, 

/ \ , ,p 2 — 2mp + e 2 , 2m 2 — e 2 p — m .„„. 

r(p)=r + p + mlog( H \ )+ tan" 1 H (24) 

& - m z Ve 2 - m 2 Ve 2 - m 2 

We choose the constant ro such that r(0) = 0. Since r is an increasing function of p this 
implicitly defines p as a function of r. It is also easy to compute from the ODE (|22|) that 

l im ^ = (3e 2 )l/3 ; lim P = h 

r— >0 t I r— >oo f 
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and result 

lim r 1 / 3 a = (e/3) 1 / 3 , lim a = 1, 

r^O r— >oo 

so that condition (iv) of Theorem |3] is clearly satisfied since 

— ► oo as r — > 0, — > 1 as r — > oo. 

ar ar 

To prove Theorem ^ it thus remains to check that the function 

1 d 2 p 2m 2e 2 + Am? 6me 2 2e 4 
p dr 2 /9 4 p 5 p 6 ' 

satisfies conditions (i-iii) of Theorem |31 In the case of condition (hi) this does not appear to be 
an easy task, because of the transcendental relation (|24|) between r and p, and the dependence 
on the two additional variables e and m. To overcome these difficulties, first we exploit the 
inherent scaling in the problem to eliminate e (effectively setting it equal to one) and then by 
forgetting the relationship between r and p and treating them as independent quantities, we 
transform the problem into one about polynomials in three variables, which we then solve by 
utilizing a few basic tools from the theory of real algebraic curves. 

Remark 1 The lower bound on the charge-to-mass ratio in Theorem ^ is not sharp. The 
actual ratio for which condition (iii) of Theorem El starts to be violated is less than 2 (It is not 
hard to see that this condition is violated for the ratios close to 1). We have not attempted to 
find the precise cut-off point, since it is not known whether condition (iii) is necessary for the 
estimate (|21|) to hold. We also note that the charge-to-mass ratios that occur in nature, such 
as in elementary particles, are in fact very large 5 . 

We perform the scaling first: Let 

r p m 2 

x := — , y := — , z := — , v := e V 

e e e 



We then have 



n 2Z 1 

a 2 = 1 + -=■. 

y v l 



Let y = 4> z {x) denote the solution to the following initial value problem 

^ = a 2 (y,z), y(0) = 0. 
The inverse function of <j> z is in fact explicit: 



= y + l °s(y 2 - tyz + 1) + 2 ! % tan^ 1 ^ 



VT^ 2 i-yz 

Note that since we are considering the super-extremal case, the parameter z ranges from to 
1. 



5 The charge e and mass m appearing in (ll'MI are in geometric rationalized units, with G = c = H = 1. The 
corresponding quantities in unrationalized electrostatic units e* and m* are related to e and m in the following 
way: m = Sjin* and e = -^e*. It follows that e/m ~ 10 21 for an electron and ~ 2 x 10 18 for a proton. 
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From the ODE it is clear that 



We have 



lim (ME)l! = 3, hi^.L (25) 

x— »0 X rr^oo x 



. 1 d 2 y 2z 2 + 4z 2 6z 2 2 . . . 2 

y dx 2 y 6 ?/ 4 y & y b y b 

From (j23|) it follows that 

2 

lim x 2 v((/> z (x),z) = --, (26) 
x^o y 

lim x 2 v{<t) z {x),z) = 0, (27) 

x— >oo 

and thus conditions (i) and (ii) of Theorem are satisfied since — | > — \ and v is clearly 
negative and increasing for y < 1/z, positive for y > 1/z, and smooth away from x = 0. 

Remark 2 Note however that since — | < |, the operator 

Q:=-d 2 x -\ + v(Mx),z) 

with domain C£°(R + ) will have non-unique self-adjoint extensions, i.e. is limit-circle at zero 
(see |2()l Appendix to X.l]). In this sense, with respect to spherically symmetric scalar waves, 
the naked singularity at x = of the super-extremal Reissner-Nordstrom solution has a £7(1) 
hair. Note also that the angular momentum contribution to the full Laplace-Beltrami operator 
of this manifold is 

a 2 {4> z {x),z) 

where £ is the spherical harmonic degree. Since this behaves like x -4 / 3 near the origin, its 
addition to Q does not change anything as far as self-adjoint extensions are concerned. It thus 
follows that with respect to scalar waves in general, this naked singularity is "infinitely hairy," 
i.e. the corresponding unitary group is certainly infinite-dimensional 6 . 

Consider now the expression in condition (iii) of Theorem [21 

x 2 — (xv((f) z (x),z)) = x 2 v(4> z (x),z) + x 3 — (4> z (x),z)a 2 ((j) z (x),z) 



where 



w(x,(j) z (x),z), 



„ ,1 ._ / 6z 1 8+28z 2 52z+32z 3 , 20+76.2 2 54z , 12\„.3 
W{X,y,Z) .- (-prH £S -b 1 ^ -y^ + ^) X 

2z 2+4z 2 1 6z 2 \ „2 



1 y y y y ' 



Let us define the function j z : IR + — > R + by 

Uy) = w(^ 1 (y),y,z). (28) 

6 This provides us with another reason for choosing the Friedrichs extension instead. 
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A long, routine computation of Taylor series shows that near y = 0, 



and 

while, for large y, 
and 



jM = l-^fy 2 + 0(y 3 ), (29) 

j'M = -^fy + o(y 2 ), (30) 

m = -^ L + 0{y). (31) 



j z (y) = -4zy~ 3 - 28zy- 4 logy + 0(y" 4 ) (32) 



j' z (y) = 12zy~ A + 112zy- 5 log y + 0(y- 5 ). (33) 



The implied constants in O-terms are locally uniform in z. Thus j z has a local maximum at 
y = 0, for < z < 1. We wish to prove that, for z < 1/2, this is in fact the global maximum 
of j z . Since it is easy to see that j z — > as y — > oo, it suffices to show that the only local 
maximum is the one at y = 0. To this end we first show that 

PROPOSITION 1 j z has no stationary points of inflection as long as z < 1/2. 

Proof: We compute 

j'ziy) = -Q^ (fe 1 (y), y' z ) + -fa 1 (y)» y» z ) a 2 ( y ^ = ^(C 1 (y) ) y^)> 



where 



312z+192z 3 

~ 1 7 

L 10+532 z 2 , 432z _ 108 \ „3 

W 8 + IT V™ ' 

24z i 64z 2 +32 120z , 18 \ _2 

3 R — — ii 7 X 

!P 2T 2T '" " 



+ 



48 \ 



and similarly we compute 



fz(y) = a2 ^ y ^ z) Q2((f> z \y),y,z), 



where 



„ /-„ „, „\ I _ 1202 i 240+1080 z 2 , -3024 z 3 -2784z 



it y° y° 

i 1360+9464 z 2 +2688z 4 , -8312 z-9856 z 3 
+ ^ + J3 

, 2200+12032 z 2 6048z , 1080 A „3 
+ ^ + ^J X 

I / 168z , -932z 2 -280 , 2072z+1216z 3 
+ [^T- + ^ + ^ 

i -3356z 2 -916 I 2688z 660 \ Jl 

H ^ 1 ^7 -yg- l X 



+ f_60z , 80+152 z 2 _ 284z , 112 \ 



X 



+ 4z _ 4 
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To prove the statement of the proposition, it is enough to show that the zero-sets of qi(x, y, z) 
and q2(x,y,z) are disjoint. The strategy is to forget about the transcendental relationship 
between x and y and treat them as independent variables, in order to be able to use results 
from the theory of real algebraic curves. 

At a stationary point of inflection both j' z and j" would vanish and hence so would the 
resultant 7 of qi,q 2 , considering both as cubic polynomials in x. Up to an integer multiple, this 
resultant is computed 8 to be 

R(qi,Q2] x) = y~ M {yz - l) 2 p(y, z), 

where 

p(y,z) = (-1536 y 9 +4608 y 7 ) z 9 

+ (3452 y 10 + 2040 y 8 - 20100 y 6 ) z 8 

+ (-3504 y 11 - 16456 y 9 - 8608 y 7 + 36120 y 5 ) z 7 

+ (1947 y 12 + 20360 y 10 + 62966 y 8 + 48272 y 6 - 33769 y 4 ) z 6 

+ (-576 y 13 - 11988 y 11 - 71800 y 9 - 153832 y 7 

- 104440 y 5 + 17900 y 3 ) z 5 

+ (72 y 14 + 3552 y 12 + 38762 y 10 + 143492 y 8 

+ 208760 y 6 + 109100 y 4 - 5530 y 2 ) z 4 (34) 
+ (-432 y 13 + 10464 y 11 - 66316 y 9 - 154672 y 7 

- 151552 y 5 - 62848 y 3 + 972 y) z 3 

+ (1152 y 12 + 15384 y 10 + 57803 y 8 + 83120 y 6 

+ 58958 y 4 + 20912 y 2 - 8l) z 2 
+ (-1440 y 11 - 10536 y 9 - 21648 y 7 - 20824 y 5 

- 11680 y 3 - 3888 y) z 

+ 720 y 10 + 2160 y 8 + 2760 y 6 + 1908 y 4 + 944 y 2 + 324. 

Note that 

qi{x, 1/z, z) = -8x 2 {z 2 - l)z 5 (2xz 3 - 2xz + l) , 
so it is zero only for x = or x = 2z \ z i • However, 

(5 z 2 — 3) z 2 

q 2 (l/(2z - 2z 3 ), 1/z, z) = ^- i— . 

1 — z z 

z = would correspond to x = y = oo, which is not of interest, and \/3/5 > \. It is thus 
enough to show that p(y, z) has no real zeros for < z < \. This is easy to establish for z = 0, 
since 

p(y, 0) = 720 y 10 + 2160 y 8 + 2760 y 6 + 1908 y 4 + 944 y 2 + 324 > (35) 
7 By definition, the resultant of two polynomials 

qi = o„n" =1 (a; - a<), qi = 6 m II£ 1 (a; - fa) 

is 

R(qi,<&;x) = a™Cnr=injLi(a l - fa) 

It can be computed from the Euclidean algorithm, or as the determinant of Sylvester's matrix or Bezout's matrix. 

8 The necessary computations here and elsewhere in the paper are performed by the computer algebra package 
PARI using exact integer arithmetic. 
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for all y. We postpone the rest of the proof of this proposition until the end of the section, and 
instead show first how this implies the desired result about V, i.e. that it satisfies condition 
(hi) of Theorem |3J 

LEMMA 5.1 The function jo has no critical points other than at y = 0. 
Proof: We have for z = 

4>q 1 {v) =y-tan~ 1 y (36) 

and 

qi (x, y, 0) = -4y- 10 x[(10y 4 + 35y 2 + 27)x 2 - 4(3y 3 + 2y 5 )x + y 6 }. 

Consider the algebraic curve qi(x, y, 0) = 0. It has a branch x = 0. The other two branches 
can be transformed into a hyperbola by changing variables to 

j_ x x 
£:=-s, i]:=-. 

y 6 y 

Let Sj denote the piece of this hyperbola that lies in the first quadrant of the (£, rf) plane, i.e. 

fi = {(e,»?) U>0,r?> 0,^,7?) =0}, 

with 

„2 | or-A . n ~ r 2 



Likewise, let X denote the part of the transcendental curve x = 4> 1 (y) in the same quadrant: 



h := 10rf + 35^7? + 27^ - 12^ - 877 + 1. (37) 



T={(£,r?) U>0,r,>0,r(C,r / )=0}, 
where r has the following expansion 

11 1 2 

- 1 + 3e-5^ + 7^ 



Recall that we have already shown by (|35[) that the function jo cannot have a stationary point 
of inflection. Therefore on each branch of the hyperbola $3 the second derivative jg must be of 
one sign. Thus one branch must correspond to minima and the other to maxima of jo- 

It is not hard to see that Sjl, the lower branch of the hyperbola fj, is well separated from 
1: Let R be the rectangle [0,1/9] x [0,6/37] in the (£77) plane. Then 

Claim 1 Sjl is contained in R while Tfl R = 0. 

Proof: The ^-intercept of Sjl is at £ = 1/9. The tangent line to at the point (1/9,0) is 
77 = — §|(£ — g). Since f)L is concave down, it lies below this tangent, and the 77-intercept of 
the tangent line is at 77 = 6/37 which shows the first part of the claim. Consider next the ODE 
satisfied by the curve x = y — tan -1 y, i.e. ^ = Ji + i ■ For y < 1 we have jj~ > ^- which upon 
integration yields 

3 

> 7T, (38) 

D 

and thus £ > |, while if y > 1 then ^| > | and thus, using the bound (|38[) at y = 1, 

/si/ x 1/1 2,1 1 

x(y) >-(y-l)+ x(l) > -(-y + -y)--> -y, 
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and thus r\ > g. Thus along X we have min(£, 77) > 1/6, which proves the second part of the 
claim. 

Let 9)u denote the upper branch of the hyperbola $). It intersects X at £ = 3 , 77 = 0. This 
point corresponds to the origin of the (x, y) plane where we know that j z has a local maximum. 
We have thus shown in the above Claim that j' cannot have any local minima. Moreover, we 
can compute the tangent lines to the two curves at (1/3,0) and we obtain that 

(dn\ 18 (drf\ 5 

Thus Sjjj is below X near (|, 0). It is easy to see that this is also the case near the 77 axis (which 
corresponds to the infinity of the (x,y) plane). Thus if X were to intersect Sju it would have 
to do so at least twice. Since a continuous function cannot have two consecutive local maxima 
without a minimum in between this is not possible and the curves Sjjj and X have no other 
common point in (the first quadrant of) the (£, rf) plane. This proves the lemma. 

■ 

The following Proposition completes the proof of Theorem ^ by showing that the quantity 
on the left in condition (iii) of Theorem El is equal to |. 

PROPOSITION 2 The function j z for < z < 1/2 has no local maximum other than at 

y = o. 

Proof: Let 

K := {(y,z)\y > 0, j' z (y) = 0, j" z {y) < 0}. 

Note that K is closed. From (|30|) we see that there is an e > such that there j' z 7^ in the 
region 0<y<e, < z < l- Similarly, from (j23J) we see that there is a Y such that j' z 7^ 
in the region Y < y < 00, < z < |. i^is therefore compact, since it is a closed subset of 
[e,F]x[0,i]. 

If K is nonempty, then there exists a point (y, z) in K with smallest z, i.e. z > z for all 
(y,z) G if. Thus j^(y) = and j"(y) < 0. By the previous lemma, z > 0. Suppose that 
7^ 0- Then by the implicit function theorem, we can solve 

9i(^ 1 (y) ) y>^) = ° 

near (y,z), to get the curve y = i?(-z). Now for 2; < z, ($(z),z) £ K so we have j z {ft{z)) > 0. 
Thus by continuity, jj(y) > 0, which leads to a contradiction. Thus we must have jf(y) = 0. 
But this is ruled out by Proposition ^ Therefore K must be empty. This proves the statement 
of Proposition El 

■ 

We conclude with the remainder of the proof of Proposition ^ First we need some defini- 
tions: For 

k,l 

a polynomial in two variables with real coefficients, let C be the curve 

C:={{y,z) £R 2 | p(y,z) = 0} 
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considered as a subset of M 2 in the usual topology. The Newton Polygon of p is defined to be 
the convex hull of the set 

N:={(k,l) eZ 2 | 

If E is an (oriented) edge of the Newton polygon with endpoints (k' E ,l' E ) and (k E ,l E ), then 
the numbers cLe, Pe and qE are defined by 



k" — k' 

d E := gcd{k E - k' E , l" E - l' E ), p E ~ 1 



d E 



QE 



I" - V 

_E_ L _E_ 

d E 



and the edge polynomial eE € TSL[t] by 



E+ i PE,l' F ;+ i 'lE 



1=0 



An edge is called outer if it maximizes some linear function ak + hi on the Newton polygon, 
where at least one of a or b is positive. 

The following compactness criteria for plane algebraic curves is proved in |22j : 

THEOREM 4 For the compactness of C it suffices that p is not divisible by y or z, and the 
edge polynomials corresponding to outer edges have no real zeros. 

The Newton polygon of the polynomial (|34[) is shown in Figure ^ It is obvious that p does 




944 1908 2760 2160 exponent of y 



Figure 1: Newton polygon of p 

not satisfy the criteria in this theorem, since some of the outer edge polynomials are of odd 
degree. To remedy this, we take a simple projective transformation of the plane and consider 
its composition with p: Let 

q (z',y'):=z' 1S p(^,^). (39) 

We need to show that q has no zeros in the region [2, oo) x (0, oo). The Newton polygon of q is 
the image of the Newton polygon of p under the linear map (k, I) >— ► (18 — k — l,k). It is shown 
in Figure [21 
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324 



powers of z' 



Figure 2: The Newton Polygon of q 

There are two outer edges, with corresponding edge polynomials 

ei(t) = 324t 10 + 994t 8 + 1908t 6 + 2760t 4 + 2160i 2 + 720, (40) 
e 2 (t) = 72t 4 - 432t 3 + 1152£ 2 - 1440t + 720. (41) 

d is clearly never zero, and the fact that e 2 has no real zeros is easily checked by Sturm's 
criterion. By the above Theorem, the zeroes of q are a compact set. If z' is maximal for this 
set then dq/dy' is zero. The resultant of q and dq/dy' with respect to y' must then be zero as 
well. This resultant is 

®1 ..J\ Jlli/J i\45/„/ i i\45/o„' T\foJ i i\f{J\3. 



where 



and 



R(q, ^;y') = z lll \z' - lf\z' + l) 45 (2z' - l){2z' + l)f{z'fg{z') 

f(z') = 15268608z /12 - 91375200z /10 + 235796896z' 8 - 336360313z' 6 
+ 278925810z' 4 - 126777721z /2 + 24542656 



g(z') = 

2238642500162400000000,2' 32 - 34444148848863120000000/ 30 
+237592851413120362800000,2' 28 - 985391370893335206960000^' 26 
+2763859141396512532788000,2' 24 - 5568163844214174878032000z /22 
+8328390177670642537641736z' 20 - 9407814473561334395291936,2' 18 
+8074051365793494550609047,2' 16 - 5249226125431353947046641,2' 14 
+2557646890465822492261299,2' 12 - 918262815118642547577717z' 10 
+238729228492213678314693z' 8 - 44988866247608231183315z /6 
+6473535589377087487753,2' 4 - 739775558130688228967,2' 2 
+ 49586421501845352448. 

/ has no real zeros at all, while all the zeros of g are contained in \z'\ < 2. Both of these facts 
are established by Sturm's criterion. It follows that the maximal value of z' is less that 2, and 
this concludes the proof of Proposition ^ 
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